Black-hole hair loss: learning about binary progenitors from ringdown signals 
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Perturbed Kerr black holes emit gravitational radiation, which (for the practical purposes of 
gravitational-wave astronomy) consists of a superposition of damped sinusoids termed quasi-normal 
modes. The frequencies and time-constants of the modes depend only on the mass and spin of the 
black hole — a consequence of the no-hair theorem. It has been proposed that a measurement of 
two or more quasi-normal modes could be used to confirm that the source is a black hole and to test 
if general relativity continues to hold in ultra-strong gravitational fields. In this paper we propose 
a practical approach to testing general relativity with quasi-normal modes. We will also argue that 
the relative amplitudes of the various quasi-normal modes encode important information about the 
origin of the perturbation that caused them. This helps in inferring the nature of the perturbation 
from an observation of the emitted quasi-normal modes. In particular, we will show that the relative 
amplitudes of the different quasi-normal modes emitted in the process of the merger of a pair of 
nonspinning black holes can be used to measure the component masses of the progenitor binary. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 

The black hole no-hair theorem states that a charged, 
stationary and axially symmetric black hole with an 
event horizon of spherical topology can be described by 
the Kerr-Newman geometry [TJ [5] , where it is character- 
ized by just three quantities - its mass, spin and electric 
charge [3HZ]. Astrophysical black holes are believed not 
to have any charge and so are described by the Kerr 
geometry [T], characterized by just their mass and spin 
angular momentum. An important consequence of the 
no-hair theorem concerns the behaviour of Kerr black 
holes when subjected to an external perturbation. There 
is strong evidence that black holes are stable against such 
perturbations [SHU] (see, however, Ref. [H]). Perturbed 
black holes regain their axisymmetric configuration by 
emitting gravitational radiation. The radiation observed 
by a detector takes the form 



h(t) 



E 



Ae mn e t/Tlmn cos (ui mn t + <j>t mn ) , (1) 



and it consists of a superposition of quasi-normal modes 
(QNMs) with characteristic mode frequencies OJZmn an d 
time-constants Tt mn . Here £ = 2, . . . , and m = —£, ...,£, 
are the spheroidal harmonic indices and n is an index 
corresponding to the overtones of each mode. The am- 
plitudes Ag mn depend on the relative orientation of the 
detector and the black hole as well as the nature of the 
perturbing agent and <fie mn are constants defining the ini- 
tial phase of the various mode^] One of the goals of this 
paper is to determine the amplitudes of the most signifi- 
cant modes excited during the merger of a pair of black 
holes. 



The mode frequencies and time-constants of a black 
hole of mass M are given by the general expressions 



Ufa 



>(j) 



M 



Ttmn = MGt mn (j) 



(2) 



1 For a recent review on black hole quasi-normal modes see Ref. 1161 



where Fi mn {j) and G£ mn (j) are functions of the dimen- 
sionless black hole spin magnitude, or Kerr parameter, 
j. All mode frequencies and time-constants then depend 
only on the mass M and the spin magnitude j of the 
black hole and no other parameter - a consequence of 
the no-hair theorem. Several authors have noted that 
this aspect of the no-hair theorem could be used to test 
if massive compact objects at galactic cores are actually 
rotating black holes described by the Kerr metric of gen- 
eral relativity [TTHT^] ; alternatively, it could be used as 
a strong field test of general relativity itself [T7j . 

The key idea behind the proposed tests is the follow- 
ing: If one can reliably decompose the observed gravi- 
tational radiation from a ringing black hole into a su- 
perposition of different modes, then the frequencies and 
time-constants of each of the modes could be used to in- 
fer the mass and spin of the black hole. If the object is 
truly a black hole, then the masses and spins obtained 
from the different modes should all be consistent within 
the measurement errors. Inconsistencies in the values 
of the masses and spins inferred from different modes 
would be an indication of the failure of general relativity 
or that the radiation was emitted from an object that 
is not a black hole. If a merging binary does not lead 
to a black hole then the inspiral phase may not result 
in a superposition of QNMs that can be characterized 
by just two parameters. Such signatures could be in- 
ferred by a model-independent analysis of the data, e.g., 
a time-frequency transform, or by assessing the posterior 
probability for alternative models. 

In this paper we estimate the relative amplitudes A^ mn 
of the various modes by fitting a superposition of QNMs 
to the radiation emitted by a merging black-hole binary 



obtained from numerical-relativity simulations. These 
simulations involve the 'simple' case of initially non- 
spinning black holes in quasi-circular orbits and differ- 
ent mass ratios of the binary, in the range 1:1 to 1:11. 
Analytical fits of the amplitudes so obtained are extrap- 
olated to mass ratios of up to 25, so as to study a variety 
of different systems. The validity of our extrapolation 
can only be confirmed by future numerical simulations of 
binary black holes with such large mass ratios. 

From the analytical fits we construct a model wave- 
form and calculate the signal-to-noise ratios (SNRs) in 
different modes during the merger of supermassive black- 
hole binaries of total mass in the range ~ 10 6 -10 8 M 
observed with the Laser Interferometer Space Antenna 
(LISArj 20 and of intermediate-mass black-hole bina- 
ries of total mass in the range ~ 100-10 3 M Q observed 
with the Einstein Telescope (ET) [2T] and advanced 
configuration of the Laser Interferometer Gravitational- 
Wave Observatory (aLIGO) 22tt25]. The response of a 
gravitational-wave detector is, of course, not separately 
sensitive to the different modes but to only their superpo- 
sition. However, it should be possible to measure the rel- 
ative strengths of the different modes by fitting a generic 
model to the observed data. We will pursue this latter 
approach in a forthcoming publication and restrict our- 
selves to a theoretical study of the relative importance of 
different modes. 

We find that over most of the parameter space 
explored, the modes with indices (£, m, n) — 
(2, 2,0), (3, 3,0), (2, 1,0) and (4,4,0) have SNRs for 
the ringdown phase larger than 500 in LISA provided 
the source is within a red-shift of z = 1 and larger than 
50 in ET provided the source is within a distance of 
lGpc. For aLIGO, the SNRs in (2,2,0) and (3,3,0) 
modes are larger than 10 in a significant region of the 
parameter space. However, other sub-dominant modes 
will not be visible in aLIGO when the source is at a 
distance of 1 Gpc or greater. In all cases black hole 
ringdown signals that result from equal-mass binaries 
can have far larger SNRs. The distance reach of LISA 
and ET to ringdown signals is large enough that one 
can expect a few events per year with quite a large 
(> 100) SNR [261 EZJ- One can, therefore, expect that 
future observations of black hole mergers will provide an 
excellent opportunity to test GR using several different 
QNMs. 

We will present a specific implementation of the test of 
the no-hair theorem and discuss a minimal and a maximal 
set of parameters that could be used to carry out such a 
test. The chief result of this paper is that the relative am- 
plitudes of the modes depend on the mass ratio q of the 



progenitor binary and that by measuring the relative am- 
plitudes, in addition to the frequency and time-constant, 
it should be possible to measure the component masses 
of the binary that led to the QNMs. 

The rest of this paper is organized as follows. Section 
|TT] is devoted to a discussion of the numerical relativity 
simulations used in this work, focusing on their accuracy, 
so as to give an idea of how reliable are our estimates of 
the relative amplitudes of different modes. The wave- 



form model used in this study is given in Section III 



stating the conventions and assumptions made in con- 
structing the model. Section |IV| constructs the ampli- 
tudes of the various modes in the ringdown signal using 
numerical simulations of Section [H] Section Hvl also deals 
with different options for identifying the ringdown phase, 
the method that was actually followed, the connection 
between the mode amplitudes and the mass ratio of the 
binary from which the black hole results and how this 
information was included in the signal model. In Sec- 
tion |VJ we discuss the detectability of the various modes 
with aLIGO, ET and LISA, and possible astrophysical 
information we can glean from such observations. Sec- 
tions |VI| and |VII| present the results from a covariance 
matrix analysis of how well we are able to measure the 
parameters of the ringing black hole and the progenitor 
binary. In Section VIII we propose a practical test of 
the no-hair theorem, making several remarks on which 
modes and parameters we could use for such a test. We 
make concluding remarks and outlook for further work 
in Section IX We use a system of units in which the 



2 Recently, the National Aeronautics an Space Administration in 
the United States opted out of the LISA mission. However, Eu- 
ropean Space Agency is pursuing an alternative that is similar in 
scope to LISA and we believe studying what science LISA could 
deliver is still very relevant. 



Newton's constant and the speed of light are both set to 
unity, c = G = 1. 



II. NUMERICAL SIMULATIONS OF MERGER 
AND RINGDOWN SIGNALS 

In this section, we shall briefly discuss how the nu- 
merical simulations were performed. We used the BAM 
code [28l [29]. The code starts with black-hole-binary 
puncture initial data (301 131] generated using a pseudo- 
spectral elliptic solver [35], and evolves them with the 
X-variant of the moving-puncture [331455] version of the 
BSSN [361 El] formulation of the 3+1 Einstein evolu- 
tion equations. We estimate initial momenta for low- 
eccentricity inspiral using the post-Newtonian methods 
outlined in [33J [33] . Spatial finite-difference derivatives 
are sixth-order accurate in the bulk 29J, Kreiss-Oliger 
dissipation terms converge at fifth order and a fourth- 
order Runge-Kutta algorithm is used for the time evolu- 
tion. Time interpolation in the Berger-Oliger-like adap- 
tive mesh refinement algorithm converges at second or- 
der accuracy. In the limit of infinite resolution the code 
is thus expected to converge with second order accuracy. 
However, in the regime of currently feasible simulations, 
the spatial finite differencing error dominates by far. Ar- 
tificial dissipation has no measurable effect on the phase 
accuracy of the waves. For well-resolved simulations we 



thus find sixth-order accuracy, as expected. 

The gravitational waves emitted by the binary are cal- 
culated from the Newman-Penrose scalar ^4, extracted 
at a distance Dl from the source. The details of our 
implementation of this procedure are given in Ref. |28j . 
Here we recall those details that are important to this 
paper. The quantity _Dl^4 is decomposed into spin- 
weighted spherical harmonics, and related to the GW 
strain as 



D L y 4 = D L (h+-ih x ) 



(3) 
(4) 



To calculate the radiated power, or luminosity, we re- 
quire h, which can be obtained by one time-integration 
of the spherical- harmonic coefficients ^^e m - In principle 
we need to fix only one constant of integration to obtain 
h from 5*4, but in practice the integration is contami- 
nated by numerical noise, see e.g., |4T)H4*2"] . As part of 
our analysis of the data, we produced both h and h, by 
(a) fixing the integration constants to ensure that h oscil- 
lates around zero, and rings down to zero, (b) removing 
low-and high-frequency noise via FFTs, and (c) remov- 
ing further spurious noise effects by subtracting low-order 
polynomial fits through the strain. 

We use results from simulations of non-spinning bina- 
ries with mass ratios q = {1, 2, 3, 4} that were previously 
presented in Refs. [39l [43] and an additional simulation 
of a q — 11 binary that was carried out as a part of this 
study. 

Detailed error analyses for the first four simulations 
were presented in Refs. [39, 43;, although those works fo- 
cused on the inspiral phase, while here we focus on the 
ringdown phase. We find that for the q = {1,2,3,4}- 
simulations, the error in the amplitude of the ringdown 
signal is dominated by the error due to wave extraction 
at a finite distance from the source. The wave extraction 
was performed at R cx = 70M for the data used in this 
paper. Waves were extracted at larger radii (80M and 
90M), but the numerical resolution was lower at these 
radii and numerical errors began to dominate the un- 
certainty of the ringdown waveform in the subdominant 
modes. We also found that the ratios of the waveform 
luminosities were remarkably robust with respect to the 
wave-extraction radius, and for these the errors are even 
lower; a similar effect was found in Ref. jUj, where the 
amplitude ratios between harmonics also played a major 
role in the study. 

The new q = 11 simulation included only two orbits 
of inspiral before merger, and was produced primarily to 
calculate the ringdown signal. The sizes of the mesh- 
refinement levels were varied to optimize both memory 
usage and numerical accuracy of the wave extraction, 
which was now performed at R cx = 100 M. Three sim- 
ulations were performed to validate the accuracy of the 
results; the resolutions at the wave extraction sphere were 
{0.533, 0.427, 0.356}/M, and the finest resolution at each 



black aim to capture the inspiral GW phase with high ac- 
curacy; our focus was on the accuracy of the ringdown. 
Here the amplitude accuracy of all the modes we con- 
sider in this paper was within 0.5%, and this uncertainty 
was dominated by the error due to extraction at a finite 
distance from the source. 

The spin of the final black hole is 0.25±0.01, and (0.3± 
0.01)% of the energy of the system is radiated during the 
last two orbits and merger and ringdown. The final black 
hole recoils by 55 ± 5 km/s. The large uncertainty in the 
recoil is due mostly to having the waveform from only 
a small number of orbits before merger, which makes 
it difficult to remove the oscillatory inspiral recoil from 
the results, as was done in Ref. [39]. Our results for this 
system are consistent with previous simulations of q = 10 
binaries 



III. 



ANTENNA RESPONSE TO A RINGDOWN 
SIGNAL 



Quasi-normal modes are transients that live for a very 
short duration in the detector band: In the case of 
intermediate-mass black holes that could be observed in 
ET, the time constant is at best about 60 ms (for a BH of 
mass 10 Mq and spin j = 0.7), while for supermassive 
black holes that could be observed in LISA the longest 
time constant is about 100 min (for a BH of mass 10 8 Mq 
and spin j = 0.7) (see Table IV). Consequently, it is not 



necessary to consider the motion of LISA or ET during 
the observation of a quasi-normal mode, at least not in 
the current evaluation of what science one might extract 
from their observation. 

Let us consider the response of an interferometric de- 
tector to a ringdown signal. We assume that the radia- 
tion is incident from a direction (9, if) with respect to, 
say, a geo-centric coordinate system. Let (e^, e^, ef ) 
be a set of orthonormal vectors representing a coordinate 
frame in which the ringdown modes take the transverse- 
traceless form; that is, the metric perturbation h%j due 
to the ringdown modes can be written in this frame as 



hij = h + el + h x e% 



(5) 



where h + and h x are the plus and cross polarizations, 
whose explicit expressions for a ringdown signal will be 
given below, and e + x are the polarization tensors given 
by 



e x = e«<g>e« + e* ®e« 



(6) 

(7) 



The response h A (t) of an interferometer, labelled by A, 
can be written as: 

h A (t) = Ff(8, <p, V) h + (t) + F*(9, <p, V) h x (t). (8) 

Here tp is the polarization angle, 

cos?/' = e e ■ ef , (9) 



and Fr x (6, <p, ip) are the antenna pattern functions of 
the detector given by, 



F A - D A p' 3 



F, 



A - T> A e j 



(10) 



where D A is the detector tensor. If e A 2 are unit vectors 
(not necessarily orthogonal to each other) along the two 
arms of an interferometer, then the detector tensor is 
given by 



D 



= e. 



e, — e 9 



(11) 



For our purposes it is most useful to express the radi- 
ation from a source in the source frame, in terms of its 
expansion in spin-weighted spherical harmonics of weight 
-2, namely ^ 2 Y lm : 



ih-, 



EE 

1=2 m=-t 



h em - 2 Y em (i ,(/)). (12) 



Here, (t, </>) refer to the co-latitude and the azimuth angle 
at which the radiation is emitted from the source; t is also 
the angle between the line-of-sight and the orbital (spin) 
angular momentum of the binary (black hole) . The com- 
plex coefficients h e,n in the expansion are referred to as 
£m modes. Explicit expressions for the first few modes in 
the post-Newtonian approximation for the inspiral phase 
of a binary's evolution can be found in Kidder 48J . It is 
useful to write the modes explicitly in terms of their real 
and imaginary parts 



■ £m 



A 



i m e 



-i<J>£ 



h*™ 



ih 



tin 



(13) 



This helps in extracting the amplitude and phase of each 
mode in terms of its plus and cross modes obtained in 
numerical simulations: 



Aim = 



^V™) 2 + (h^y, $ lm = tan" 



,£m 



hf 1 



(14) 

Noting that - 2 Y lm (i, 4>) = - 2 Y em (i, 0) e"™^, we can re- 
arrange the sums in Eq. ( 12 ) using Eq. ( 14 1 to get [40] 

h+= J2 A ^ Y+ l cos($, m - m<j>), (15) 

£,m>0 

^x = - E A ^ F x m sin($, m - m4>), (16) 

where we have dropped the "memory-effect" m = 
terms, for which the amplitude is low (see, e.g., a recent 
numerical study [IH]). Note that while these modes are 
non-oscillatory during inspiral, they do exhibit ringdown, 
which has been studied in some detail with numerical 
codes in axial symmetry, where they are the only non- 
zero modes (see, e.g., [SO])- In the above expressions, the 
angular functions Y^\ (t) are the following combinations 
of the spin- weighted spherical harmonics: 

y* m (0 = ~ 2 Y em (L, o) + (-i)' - 2 Y e - m (i, o), 

Y* m ( t ) = - 2 Y im (i,Q) - {-if _ 2 Y e - m (i,0). (17) 



For the inspiral phase of a binary when the two com- 
pact bodies are widely separated, post-Newtonian (PN) 
approximation gives the amplitudes Ag m {t) and phases 
3>e m (t) as expansions in v/c, where v is the velocity of the 
bodies (see, e.g., [H]). Numerical relativity simulations 
can be used to extract them when the PN approximation 
breaks down. In the case of perturbed black holes, which 
a binary will result in, black hole perturbation theory 
predicts that the modes are damped sinusoids with their 
amplitudes and phases given by: 



a im M p _ t/Ten 



Mm 



D, 



$em{t) = UJtmt, 



(18) 



where M is the mass of the black hole and Dl is its lu- 
minosity distance from Earth. Time-constants Tg m and 
frequencies Lu# m can be computed from black hole pertur- 
bation theory (see, e.g., Ref. [18] for a recent comprehen- 
sive listing of frequencies and time-constants). However, 
the amplitudes ai m depend on the nature of the pertur- 
bation and are not analytically accessible in the case of 
black holes that form from the coalescence of a binary. 
We shall "measure" them later in this paper using results 
of our numerical simulations. 

Using the above equations, the output of the numeri- 
cal simulations for the plus and cross polarizations cor- 
responds to the following expressions: 



h+(t) = ]T a jy J - Yl m (i) e-'/ r *» cos (w tm t - mfl , 

i,m>0 L 

hx (t) = J2 ^TT- Y * m W e ' t,Tlm sin ^ lmt ~ m ® ■ 



£,m>0 



D h 



(19) 



We have dropped the overtone index n from all the 
relevant quantities (amplitudes, frequencies and time- 
constants) of quasi-normal modes, as we are assuming 
that higher (i.e., n > 0) overtones, quickly become negli- 
gible in amplitude, compared to the fundamental n = 
overtone. Only the fundamental n = overtone of the 
various modes is considered in this paper. 

The amplitudes Ag m of the various modes depend on 
the nature of the perturbation. For ringdowns resulting 
from the merger of two non-spinning black holes, Ag m 
depend on the mass M of the final black hole and mass 
ratio q = m\jra 2 (q > 1) of the progenitor binary. 

In the next section we will estimate the amplitude 
of dominant modes by fitting the late time signal from 
a numerical relativity simulation to a superposition of 
ringdown modes. For binaries with non-spinning black 
holes considered in this paper we find that modes with 
(£,m) = (2, 2), (2,1), (3, 3) and (4,4) are the ones that 
are excited with amplitudes large enough to be interest- 

Ira 



ing. For these modes, the angular functions Y^ x (l) are 



given by 



yf(0 



Y?(l) = 
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cos t sin l, 



21 (1 + cos 2 t) 



costsmt, 



stant phases of each quasi-normal mode [J The effective 
amplitudes Bi m are proportional to the intrinsic ampli- 
tudes ai m of the modes and vary inversely with the lu- 
minosity distance. Their magnitude also depends on the 
various angles (6, tp, ip, t) describing the position of the 

(20) source on the sky and its orientation relative to the de- 
tector through the antenna pattern functions F + and F x 
and spherical harmonic functions Yf m and Y^ m . The 
constant phases j£ m also depend on the angles and the 
fiducial azimuth angle cf>. 

The above form of the response is more helpful in un- 

(21) derstanding which, or which combination, of the param- 
eters can be measured and how many detectors are re- 
quired in solving the inverse problem, namely to fully 
reconstruct the incident gravitational wave and the pa- 
rameters of the source that emitted the radiation. We 
shall use the above form of the waveform to compute the 

zoo") signal-to-noise ratios and the covariance matrix. 



Yt%) 
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Tott 
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sin 2 i, 



cost sin /,. 



(23) 



In the case of binaries comprising of black holes with 
generic spins, the relative amplitudes of the various 
modes will also depend on the magnitude and direction of 
the spin vectors of the progenitor black holes. A detailed 
study of the dependence of the relative amplitudes of the 
various modes on the initial spin configurations and mass 
ratio of the progenitor binary is necessary to assess how 
accurately one might be able to use quasi-normal modes 
to measure a progenitor binary's parameters. For this, a 
more exhaustive set of simulations covering the full pa- 
rameter space of binary black holes is required and will 
be taken up in the future. 

We conclude this section by noting that using the ex- 
pression for the two polarizations in Eq. ( 19 1 the detector 
response given in Eq. (T8J) can be written as: 



h: 



e,m>o 



-t/n m cos ( Utm t + 7ftn ) 



(24) 



where the superscript A is an index denoting the detector 
in question (which is relevant when we have a network 
of detectors), Bg m (ji m ) is the following combination of 
the amplitudes Ag m (respectively, phases mcf>), antenna 
pattern functions F+ and F x and the inclination angle 



Bfrn, = 



llv 



^L^ {F Ay im) - + {F A Y e. m) ^ (25) 

(26) 
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b( m + m <p + tan 



J? A y£m 
pi 'yttn 



Note that, for the sake of clarity, we have dropped the 
index A on Bg m and "fi m - Here, (j)£ m are arbitrary con- 



IV. AMPLITUDES OF MODES EXCITED 

DURING THE RINGDOWN PHASE OF A 

BLACK HOLE BINARY 

In this section we will use numerical simulations to 
evaluate the amplitude of the various modes excited as 
a function of the mass ratio q. We will examine how the 
amplitude of the dominant 22-mode varies as a function 
of q. Of particular interest would be the growth of 21-, 
33- and 44-modes relative to the 22-mode as the binary 
system becomes more asymmetric. We will provide sim- 
ple analytical fits to the amplitudes of all the different 
modes and discuss how the inclination angle of the black 
hole's final spin will affect the amplitude of the ringdown 
signal. Our analysis is complementary to previous stud- 
ies of the mode structure of unequal-mass nonspinning 
binaries, for example [301 IH3| • 



A. Evolution of the luminosity 

An important question that arises in the study of 
QNMs excited during the merger of a black hole binary 
is the determination of the most dominant modes in the 
infinite mode sum in Eq. (fTl) for this particular kind of 
perturber. This analysis is necessary to construct a good 
model of the waveform to use in the analysis. 

Figure [l] plots the luminosity in gravitational waves, 



L 



tin 



= Dl 



£m\2 



(hf 1 ) 2 + (hf 1 ) 



in the first four most dom- 



inant modes as a function of the dimensionless time t/M, 
where M is the initial total mass of the binary. The lu- 
minosities are plotted for four values of the mass ratio 



3 Specifically, the £ = m modes have a nearly consistent rotational 
phasing, while the £ j^ m modes seem to have somewhat distinct 
associated dynamics, with differentiated amplitude and phasing 
relationships during the merger process. |52| 
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FIG. 1. This plot shows the relative luminosities, or radiated power, in modes (2,2), (3,3), (4,4), and (2, 1). L; m represent the 
luminosities (in units c = G = 1 in which luminosity is dimensionless) and ri m denote the ratios r; m = Li m /L22- The different 
panels correspond to systems with different mass ratios as indicated in the panel. Note that as the mass ratio increases, the 
luminosity in each mode decreases but the amplitudes of all higher-order modes relative to the (2, 2)-mode increase. We have 
omitted — both in the figure and in this work — the next most dominant modes, (5, 5), (3, 2), (4, 3), (6, 6) and (5, 4) as they 
are generally less than one percent as luminous as the (2, 2) mode, (see, however, Pan et al [51]). 



q = 2,3,4 and 11. Wc have left out the plot correspond- 
ing to the equal mass case q — 1, as in this case the modes 
with odd values of I or m are not excited and hence not 
as interesting as when the masses are unequal. The lu- 
minosity peaks just before the two black holes collide 
but different modes peak at different times. The 21- and 
33-modes peak after the 22-mode reaches its maximum. 
However, the 44-mode shows the opposite behaviour. For 
a more thorough investigation on the different multipolar 
contributions to the total radiated energy, see Ref. [52"] . 

Although the 33-mode is absent when q = 1, it is al- 
ready more dominant than the 44-mode when q = 2 and 
remains the most dominant after the 22-mode, through- 
out the inspiral and merger phase and for all mass ratios 
(except, of course, when q = 1). The 44-mode remains 
more dominant than the 21-mode for the most part, but 
the trend reverses after merger. This is because the 21- 
mode reaches its peak luminosity a little after the 44- 
mode. For q — 2, when the 22-mode reaches its peak, 
the luminosity in the 33-mode is an order-of-magnitude 



smaller than the dominant 22-mode; luminosities in 21- 
and 44-modes are 50 times smaller than the 22-mode. 



In addition to the luminosities, we have also plotted 
their ratios rg rn — Le„J L22 with respect to the 22-mode. 
It is clear that for more symmetric systems (i.e. q <; 1) 
higher modes are hardly excited. For instance, when 
q = 2 the luminosity of the 21- and 44-modes remains 
below ~ few percent of the 22-mode, while the 33-mode 
is always less than ~ 5% of the 22-mode. As the mass 
ratio increases the higher order modes are excited with 
greater amplitudes. The different mode amplitudes be- 
come comparable to the 22-mode and to one another as 
the mass ratio increases. In the next two sub-sections we 
will give a more quantitative evaluation of the relative 
mode amplitudes in the ringdown part of the signal. 




FIG. 2. Evolution of the first few dimensionless mode frequencies fi m = Mu)i m as a function of the dimensionless time t/M, 
for different values of the mass ratio q of the progenitor binary. Also shown in arbitrary units is the luminosity in the 22 mode. 
All mode frequencies, especially /22 and /33, stop evolving and stabilise soon after the binary merges to form a single black 
hole. The waveform is assumed to contain a superposition of only quasi-normal modes a duration 10M after the luminosity in 
22 mode reaches its peak. 



B. Identifying the ringdown phase 

An important task of our study is the identification 
of the point when the signal is purely a superposition 
of various quasi-normal modes. This is necessary in or- 
der that the proposed tests of general relativity are not 
corrupted due to the presence of extraneous signals. By 
assuming that the ringdown phase occurs sooner than it 
actually does, we are in danger of corrupting the wave- 
form. Equally, by identifying the ringdown phase to be 
much later than it actually is, we will significantly weaken 
the tests since the signal amplitude falls off exponentially 
from the beginning of the ringdown phase. A proper 
identification of the beginning of the ringdown phase is 
needed to correctly extract the amplitude of the quasi- 
normal modes and to compute the signal-to-noise ratio 
and other quantities. 

To this end, we shall use the evolution of the frequency 
of the various modes hi m . As a binary evolves, the fre- 
quency of each mode hi m increases, the rate of increase it- 
self being greater as the two black holes get closer. When 
the two black holes merge, a common horizon forms and 
the frequency of each mode stabilizes, finally reaching the 
quasi-normal mode value as predicted by black hole per- 



turbation theory. We shall identify the beginning of the 
ringdown phase to be (approximately) the epoch when 
the frequency of the various modes begin to stabilize. 

We can compute the frequency of each mode from 
the evolution of its phase given by the second of the 
equations in Eq. (14). Once the phase is known it is 



straightforward to write down the (dimensionless) fre- 
quency f im = Mut m = d$t m (t)/d(t/M). The ringdown 
phase can be assumed to begin when fi m computed from 
our numerical simulations are close to those obtained 
from black hole perturbation theory. We will first take 
a look at the predictions from black hole perturbation 
theory and then compare those predictions to the results 
obtained from our numerical simulations and plotted in 
Fig.[2] 

There has been a lot of work on the computation of 
the frequencies and time-constants of various modes of a 
perturbed Kerr black hole. Berti et al [TH] have found 
simple fits, as a function of the spin parameter j, to the 
dimensionless mode frequencies/] fi m — M ug m and qual- 



Note that we only consider here, the least damped n = overtone 



TABLE I. Dimensionless frequencies fi m for various modes 
for different values of the black hole spin j that results from 
the merger of a binary of mass ratio q. We assume that the 
ringdown phase begins when the frequencies of various modes 
stop increasing and stabilize to a constant value. For each 
mode frequency, the column labelled "Fit" gives the values 
computed using Eqs. ( 2T l-( 30 1 and that labelled "NR" shows 



values at the "beginning" of the ringdown mode, a duration 
10M after the gravitational wave luminosity in the 22 mode 
reaches its peak (dashed vertical line in Fig.[2|. 



Q 


3 


/22 

Fit NR 


Fit NR 


/33 

Fit NR 


/44 

Fit NR 


1 


0.69 


0.53 0.51 


0.46 - 


0.84 - 


1.14 1.08 


2 


0.62 


0.50 0.49 


0.44 0.42 


0.80 0.78 


1.09 1.05 


3 


0.54 


0.48 0.47 


0.43 0.41 


0.76 0.74 


1.03 1.01 


4 


0.47 


0.46 0.45 


0.42 0.43 


0.73 0.72 


0.99 0.97 


11 


0.25 


0.41 0.39 


0.39 0.41 


0.66 0.64 


0.89 0.85 



ity factors 2 Qe m — i~i m uj£ rn . The fitting functions for the 
22, 21, 33 and 44 modes are given by [18] 



J22 = 1-5251 - 1.1568(1 - j) 01292 , 
Q 22 = 0.7000 + 1.4187(1 - j)- 0A990 , 



(27) 



f 21 = 0.6000 - 0.2339(1 - j) 04175 , 
Q 21 = -0.3000 + 2.3561(1 - j) -02277 , (28) 



/33 = 1.8956- 1.3043(1 -j 

Q 33 = 0.9000 + 2.3430(1 - j 



jO.1818 
1-0.4810 



744 

Qi4 



2.3000- 1.5056(1 -j) 02244 , 
1.1929 + 3.1191(1 ^-°- 4825 



•3) 



(29) 



(30) 



These fits are quite robust and they differ from the actual 
values obtained for the frequencies and quality factors by 
no more than 3% [IS]. Table [I] lists frequencies /22, hi, 
/33 and /44, for several different values of the spin param- 
eter j. The chosen values of j correspond to the final spins 
of black holes that result in our numerical simulations of 
binaries with different mass ratios q. Values in columns 
labelled "Fit" are those obtained using Eqs. ( 27 1-( 30 ) and 



those in columns labelled "NR" are those obtained from 
our numerical relativity simulations as follows. 

Figure [2] plots the frequencies fe m for the 22, 21, 33, 
and 44 modes as (black) dot-dashed, (blue) dot-dot- 
dashed, (green) dashed and (red) dotted curves. As 
expected, the frequency of each mode increases, quite 
rapidly towards the end, but stabilizes to a constant value 

the quasi-normal mode frequency of the final black hole. 
The plots also show the (arbitrarily scaled) luminosity in 



for each mode and have therefore dropped the overtone index 
from mode frequencies, quality factors and time-constants. 



the 22-mode as a (black) solid curve. The epoch at which 
the luminosity reaches its peak (indicated by a solid ver- 
tical line) has been set to be t/M — 0. We see that the 
onset of the ringdown phase occurs significantly after the 
luminosity reaches its peak. For simplicity, we have cho- 
sen the beginning of the ringdown phase to be a duration 
t = 10 M after the system reaches its peak luminosity, in- 
dicated by a dashed vertical line in Fig. [2] In reality for 
the I — m modes, it is 4M-5M earlier than for the (2, 1) 
mode. The luminosity curves, Fig.[l] exhibit a similar 
behaviour. See also 52 . 

Most modes seem to stabilize at the onset of the ring- 
down phase. For q ;$ 4 the frequency of 22 and 33 modes 
stabilize, but it is less so with 21 and 44 modes. In our 
q = 11 simulations all mode frequencies seem to oscillate 
around a mean value, 22 and 21 more than 33 and 44. 
From Fig. [I] we see that after reaching the peak lumi- 
nosity the amplitudes of 21 and 33 modes relative to the 
22 mode (thick lines) are constant. This justifies why we 
might fit a g-dependent function to the relative ampli- 
tudes of various modes (see Sec. |IVC and Fig. [3]) that is 
valid throughout the ringdown phase. The same cannot 
be said about the 44 mode. 

Under columns labelled "NR" , Table [I] gives frequen- 
cies / 22 , /2I) J33 and / 4 4 at the onset of the ringdown 
mode (i.e., an epoch t = 10M after the luminosity of the 
22 mode reaches its peak) obtained from our numerical 
relativity simulations. Modes with odd £ or m are not 
excited when a binary comprising a pair of equal mass 
black holes merges, which is the reason why these entries 
are missing from the Table. The mode frequencies at 
10M after peak luminosity agree with the fits obtained 
from black hole perturbation theory to within 5%. Hence 
we believe that our method is quite robust in identifying 
the ringdown phase. 



C. Relative amplitudes in the ringdown phase 

Our goal here is to estimate the amplitudes of the var- 
ious modes in the ringdown phase. The transition from 
the inspiral to the ringdown phase is very smooth and it 
is not easy to pick a unique instant after which the tran- 
sition occurs. Amplitudes of the various modes in the 
ringdown phase are given in Table III] at three different 
epochs: 

• at the epoch when the luminosity of the 22 mode 
reaches its peak, 

• a duration 10M after the 22 mode reaches its peak 
luminosity, and 

• a duration 15M after the 22 mode reaches its peak 
luminosity. 

The Table lists the absolute amplitude a 22 of the 22 mode 
and relative amplitudes a^ m /a 2 2 of the rest of the modes. 
These amplitudes are plotted in Fig. [3] 



TABLE II. For different mass ratios (column 1), we list the final spin j of the black hole (column 2), amplitudes of different 
modes at three different epochs: (i) when the 22 mode reaches its peak luminosity (columns 3-6), an epoch 10M (columns 
7-10) and 15M after the 22 mode reaches its peak luminosity. We list the absolute value of the amplitude for the 22 mode 
and the ratio of amplitudes of the rest of the modes to the 22 mode. The beginning of ringdown, taken as the point when the 
instantaneous frequency of each mode begins to stabilize, is typically found to be ~ 1-2 cycles after the peak luminosity. For 
concreteness we take the beginning of the ringdown mode to be 10M after the peak luminosity. 



<? 


3 


At peak luminosity 


At 10A/ after peak 


At 15M after peak 






«22 


033/022 044/022 


"21/022 


«22 033/022 Q 44/Q22 


021/022 


Q22 033/022 044/022 


021/022 


1 


0.69 


0.365 


0.000 0.052 


0.000 


0.217 0.000 0.043 


0.000 


0.152 0.000 0.038 


0.000 


2 


0.62 


0.321 


0.149 0.050 


0.114 


0.194 0.161 0.030 


0.121 


0.132 0.156 0.020 


0.135 


3 


0.54 


0.266 


0.216 0.070 


0.178 


0.158 0.247 0.052 


0.195 


0.108 0.244 0.044 


0.212 


4 


0.47 


0.225 


0.259 0.087 


0.203 


0.140 0.267 0.069 


0.234 


0.095 0.264 0.057 


0.262 


11 


0.25 


0.100 


0.349 0.156 


0.312 


0.063 0.377 0.154 


0.407 


0.048 0.347 0.137 


0.478 



First let us note that although, as expected, the ab- 
solute amplitude of the 22 mode depends on the epoch 
that we identify as the start of the ringdown phase, the 
amplitudes of the sub-dominant modes relative to the 22 
mode are not too sensitive to that identification. This 
is especially true for the 33 and 44 modes whose peak 
luminosity is at the same epoch as that of the 22 mode 
[thin (green) dashed and (red) dotted curves in Fig. [I] , 
but less so for the 21 mode [thick (blue) dot-dot-dashed 
curve in Fig. [Tl bottom right panel of Fig. pj], whose peak 
luminosity occurs significantly after that of the 22 mode. 
As mentioned before, for concreteness we shall take the 
starting point of the ringdown mode as an epoch lOAf af- 
ter the 22 mode reaches its peak luminosity (the dashed 
vertical line in Fig. [2]). All discussions in the remainder of 
this paper are based on this identification. A duration of 
10M corresponds to between 1 and 2 gravitational-wave 
cycles of the merger signal. 

As the mass ratio of the progenitor binary increases, 
the amplitude of the 22 mode rapidly decreases but the 
sub-dominant modes approach each other in power (see 
Fig. ni and their amplitudes increase (cf. Table In] and 
Fig. p|. For a mass ratio of q = 4, the 33 and 21 modes 
have amplitudes 1/4 that of 22 while at q = 11 they are 
40% of the 22 mode. Of course, the overall luminosity 
decreases as the mass ratio increases and one expects no 
radiation in the limit q — > 00. Indeed, the emitted energy 
during the 'merger' phase goes roughly as the square of 
the symmetric mass ratio v — mim2/M 2 of the progeni- 
tor binary [4"0] . 

Let us note that the values in TableHJdo n °t all refer to 
the same final black hole spin. All our black holes are ini- 
tially non-spinning and the final spin is simply the resid- 
ual angular momentum of the progenitor binary. The 
final spin, therefore, depends on the mass ratio, and is 
greatest when the two black holes are of the same mass. 
In principle, it should be possible, but in practice very 
difficult, to produce numerical data for different mass ra- 
tios all with the same final spin. To do so we require 
an accurate mapping between the mass ratio and ini- 
tial black-hole spins, and the spin of the final black hole. 



For example, configurations that lead to a nonspinning 
Schwarzschild black hole are suggested in Ref. [S3] and 
more generic cases are considered in Refs. [5"5Tl58| . How- 
ever, fine-tuning the spin of the final black hole requires 
that the component black holes are also spinning, and in 
this paper we consider only binaries with non-spinning 
components. 



D. Fitting functions for relative amplitudes 

For the purpose of computing the signal-to-noise ratio 
and the covariance matrix it is convenient to have analyt- 
ical expressions for the relative amplitudes of the various 
modes. Our fits are meant to capture the dependence of 
the amplitudes on the mass ratio in the range we have 
considered in this paper. They are not meant to explore 
the complex dynamics of a ringing black hole and are not 
necessarily valid outside the region we have explored. 

We have seen that the amplitudes inferred from the 
simulations depend on how we identify the beginning of 
the ringdown phase. As mentioned before, in our calcu- 
lations we have assumed that the ringdown phase begins 
an epoch 10M after the 22 mode reaches its peak lumi- 
nosity. We find that the amplitude of the 22 mode at this 
epoch as a function of the mass ratio is well fitted by 

a 22 (q) = 0.25 e-*/ 7 - 5 . (31) 

The solid curve in the top left panel in Fig. [3] shows that 
this is a good fit to the data points (filled squares). Dot- 
dot-dashed and dotted lines in the same panel show the 
fits that describe 0:22 (<?) at the peak of the luminosity and 
15M after the peak, respectively. In all cases, the data 
is pretty well approximated by an exponentially falling 
function. 

The relative amplitudes of the sub-dominant 21, 33 
and 44 modes are well fitted by the following functions 

a 2 i(<7)=0.13a 2 2(<z)(<Z-l) 1/2 , (32) 

a 33 (q)= 0.l8a 22 (q) (q - 1) 1/3 , (33) 

a A4 {q)= 0.024 a 22 {q)q 3/4 . (34) 
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is given by (see, e.g., [55] ) 



FIG. 3. This plot shows the amplitudes as a function of the 
mass ratio for different modes at the peak of the luminosity of 
the 22 mode (circles), an epoch 10M and 15M after the peak 
(respectively, squares and triangles). We have plotted the 
absolute amplitude Q22 of the 22 mode and ratio of the sub- 
dominant mode amplitudes ctt m /a.22 relative to 22 (cf. Table 
IHb. The solid lines are the best fits [cf. Eqs. <[3l])-<[34J)] to the 
amplitudes at 10M after the peak luminosity. 



They are plotted as solid curves in the relevant panels 
of Fig. [3j The fits are motivated by the fact that when 
q = 1 only modes with even values of I and m are ex- 
cited, while those with odd I or in are absent. Note that 
the relative amplitudes grow as a binary becomes more 
asymmetric and so higher order modes should be more 
easily detectable if the mass ratio is large. In absolute 
terms, of course, all modes are exponentially damped as 
a function of q. 

In this work we estimated the relative amplitudes of 
different modes at a single epoch. It might be more reli- 
able to estimate their average over a small duration, for 
instance between 10M-15M. This could diminish any 
'numerical noise' present and extract more accurate fits 
but we did not explore this alternative approach in this 
work. 



V. VISIBILITY OF RINGDOWN MODES 

In this section we will study the signal-to-noise ratio 
(SNR) obtained by various detectors for the ringdown 
phase of the coalescence of a black hole binary. We will 
begin by defining the matched filter SNR, followed by the 
noise power spectral densities of aLIGO, ET and LISA 
and the choice of signal parameters used in the study. 
We will then discuss the visibility of the different modes, 
focusing on the distance reach of the various detectors. 



A. Matched filter SNR 

The matched filter SNR p obtained while searching for 
a signal of known shape buried in Gaussian background 



^ 1 / ™^df 
" 4 /o W d/ 



(35) 



where Sh(f) is the detector noise power spectral density 
and H(f) is the Fourier transform of the signal assumed 
(in this work) to be a superposition of the 22, 21 and 
33 modes. In the time-domain our waveform is given by 
Eq. (24]) where the sum is over (£, m) = (2, 2), (2, 1), (3, 3) 
and we have ignored all higher order modes including the 
(4, 4) mode. The coefficients ag m required to compute 
the waveform are assumed to be as in Eqs. ( 31 )-( 33 ). 



The ringdown signal is a superposition of different 
modes and its visibility depends not only on the relative 
amplitudes but also on the relative phases of the different 
modes. To this end, it is useful to define the SNR of a 
mode as 



2 

Ptn 



Sh{f) 



d/, 



(36) 



with the caveat that the total SNR p 2 is not the quadra- 
ture sum of p 2 m since there are also interference terms 
that can be negative. 



B. Sensitivity curves 

In our study, we will consider the performance of three 
detectors: the aLIGO ET and LISA. A fit to the aLIGO 
noise spectral density tuned to detect binary neutron 
stars iqj 



S h (f) = 10 



-49 



123.35 



10 16-4(/-7.9) 2 +008 ^4.69 

1 -0.23 x 2 + 0.0764 x 4 



1 + O.Ux 2 



Hz" 1 (37) 



where x = //215 Hz. In the case of ET we consider the 
sensitivity curve designated ET-B [60] whose noise power 
spectral density is given by Sh(f) = 10~ 50 h n (f) 2 Hz -1 



h n (f) = 2.39 x 10 
+ 1.76 x" 



-27^-15.64 



0.349 x 



-2.145 



0.409 x 



1.10 



(38) 



where x = //100 Hz. We take LISA noise spectral density 
to be the one that was used by the LISA Parameter Esti- 
mation Taskforce in Ref. [61] , which also corresponds to 
the noise curve from the second round of the Mock LISA 
Data Challenge [51 [53]. 



The fit was provided by C. Capano, Syracuse University. 
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FIG. 4. The signal-to-noise ratio integrand for LISA for a quasi-normal mode signal that is composed of 22, 21 and 33 modes 
— the three most dominant ones. The source is assumed to be at a red-shift of z — 1 and the various angles are as in Tabic [ITT] 
The left panel corresponds to a black hole of mass M — 5 x 10 6 Mq and the right panel to a black hole of mass M — 10 7 Mq. 
In both cases the mass ratio of the progenitor binary is taken to be q = 10. 



C. Choice of various parameters 



The SNR depends on a number of source parameters 
as well as the location of the source on the sky. We 
have limited our investigations to studying the SNR and 
covariance matrix as a function of the black hole's (ob- 
served) mass M and the mass ratio q of the progenitor 
binary, for fixed values of the distance to the black hole 
and various angles. In the case of LISA the black hole 
is assumed to be at a red-shift of z = 1 which corre- 
sponds (in our cosmological model) to a luminosity dis- 
tance of Dl — 6.73 Gpc. In the case of aLIGO and ET 
we set £>l = 1 Gpc. In all cases the angles are fixed to 
be 6 = ip = l = ip — 7r/3. For a statistical analysis of the 
effect of these angular parameters on the detectability of 
a ringdown signal see |64) . 

The black hole mass is varied over the range 
[100, 1O 3 ]M in the case of aLIGO, [10, 10 3 ] M in the 
case of ET and [3 x 10 6 , 10 8 ] M© in the case of LISA. 
These choices are dictated by the frequency sensitivity of 
the instruments which further dictates the range of black 
hole masses whose ringdown radiation they are most sen- 
sitive to. 



TABLE III. This Table lists the values of the various param- 
eters used in our study. In all cases the angles are all set to 

8 — tp — Tp = l = (p = 7r/3. 



Our choice of parameters is summarized in Table |IH| 
We reiterate that our masses are observed masses which 
means that the intrinsic mass of the black hole is smaller 
by a factor 1 + z ~ 1.2, for aLIGO and ET and by a factor 
of 1 + z = 2 in the case of LISA. Although the signal 
visibility simply scales as the inverse of the distance, the 
fact that the mass is blue-shifted means that we cannot 
easily scale our results to another (say, a greater) distance 
for the same intrinsic masses. Such a scaling will be valid 
if at the same time the intrinsic masses are also scaled 
up/down by the appropriate red-shift factor. 



Detector D^/Gpc 


M/Mq 


q 


aLIGO 1.00 


[100, 10 3 ] 


2-10 


ET 1.00 


[10, 10 3 ] 


2-10 


LISA 6.73 


[3 x 10 6 , 10 8 


2-25 



D. Visibility of different modes 

It is instructive to plot the SNR integrand dp 2 /df = 
\H(f)\ 2 /Sh(f) as it depicts how the different modes be- 
come important for systems with different masses. Fig. HI 
plots this quantity for two systems as seen in LISA. 
The left panel corresponds to a black hole of total mass 
M = 5 x 10 6 M© and the right panel to M = 1O 7 M . 
The mass ratio q is q = 10 in both cases and the angles 
are as in Table Mil 

The various mode frequencies F(. m — ujf m /(2ir) of the 
two systems are F22 — 2.74 mHz, F21 — 2.54 mHz and 
F33 ~ 4.26 mHz for the lighter black hole and F22 — 
1.37 mHz, F 2 i ~ 1.27 mHz and F 33 ~ 2.13 mHz for the 
heavier black hole. Let us first note that the 22 mode of 
the lighter black hole and 33 mode of the heavier black 
hole are close to the region where LISA has best sensitiv- 
ity. This will be relevant in the discussion that follows. 

The intrinsic amplitudes of the 21 and 33 modes are 
a little more than a third of the 22 mode for q = 10. 
However, since the SNR integrand depends on the signal 
power weighted down by the noise power, for a given 
black hole mass the SNR integrand could be as large 
as, or even dominated by, modes different from the 22 
mode. This does not happen for the 21 mode since the 
frequencies of the 22 and 21 modes are very close to each 
other and so the 21 mode is always far smaller than the 
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FIG. 5. Signal-to-noise ratio (SNR) in Advanced LIGO (top set of four panels) and Einstein Telescope (bottom set of four 
panels) as a function of the black hole's mass M and progenitor binary's mass ratio q for different modes. Most of the 
contribution to the SNR comes from the 22 mode but other modes too have significant contributions, 33 being more important 
than 21. The source is assumed to be at a distance of 1 Gpc and various angles are as in Table [Till 
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22 mode. For a black hole mass of 10 7 M the 33 mode 
is as strong as the 22 mode and for masses even larger, 
the 33 mode overwhelms the 22 mode. The total SNR for 
the 5 x 10 6 M Q black hole is p = 1670, with the different 
modes contributing P22 = 1500, P21 — 625, and P33 = 
950. The SNR is clearly dominated by the 22 mode. 

In the case of the heavier black hole, the total SNR is 
p = 2520, with the different modes contributing p 2 2 — 
1940, P21 = 920, P33 = 1860. In this case, the 33 mode is 
as strong as the 22 mode but the 21 mode, as expected, 
is sub-dominant. 



E. Exploring black hole demographics with ET and 

LISA 



Figures [5] and [6] plot the SNR in the ringdown signal 
(plots titled "SNR in all modes" ) and contribution from 
the 22, 21 and 33 modes (plots titled accordingly) as 
a function of the black hole mass M and mass ratio q 
of the progenitor binary for aLIGO, ET and LISA; M 
and q are varied over the range as in Table |Pll| Most of 
the contribution to the SNR comes from the 22 mode 
followed by 33 and 21. Let us recall that SNRs from 
different modes do not add in quadrature. 

In the case of aLIGO, the 22 and 33 modes will be vis- 
ible in a significant fraction of the parameter space ex- 
plored provided the source is within a distance of 1 Gpc. 
The 21 mode will not be visible in aLIGO at this distance 
except perhaps for the heaviest systems explored. 

In the case of ET, assuming a SNR threshold of 10 for 
detection, the signal is visible to a red-shift of z ~ 0.8 
in most of the parameter space explored. Black holes of 
total mass M > 400M Q that form from the coalescence 
of binaries whose mass ratio is less than 4 will be visible 
at red-shifts z ~ 2-3.5 

In the case of LISA, ringdowns produce a very large 
SNR. Even assuming an SNR threshold of 40, LISA 
should see the formation of supermassive black holes in 
the range [10 6 , 10 8 ] M Q up to a red-shift of at least z ~ 6 
but if the progenitor black holes have mass ratio q < 10 
they should be visible from the earliest moments of their 
formation in the Universe. 

Our results unambiguously demonstrate that ET and 
LISA can together probe black hole demographics, ET 
exploring the lower end of the mass spectrum of seed 
black holes and LISA the higher end of that spectrum. 
The two detectors together cover a large mass range from 
~ 10 2 M Q , all the way to ~ 10 8 M & . (Although out of the 
range of masses explored, note that ET could observe 
heavier black holes of mass 10 4 Mq and LISA could ex- 
plore lighter black holes of mass 10 5 M Q .) The distance 
reach will be different depending on the mass ratio of 
the progenitor binary and the total mass of the black 
hole. Even so, LISA and ET will make it possible to ex- 
plore the formation of black holes and trace their merger 
histories and possibly help understand the role of black 
holes as seeds of galaxies and large scale structure in the 



Universe. 



VI. WHAT CAN A RINGDOWN SIGNAL 
MEASURE? 



By measuring the ringdown signal and resolving it into 
different modes, we should be able to learn a great deal 
about the merger dynamics and test general relativity. 
For instance, by determining the total mass of the binary 
from the inspiral phase and comparing it to the mass of 
the black hole obtained from the ringdown phase we can 
measure, quite precisely, how much mass is converted into 
radiation in the process of merger. LISA can typically 
measure the total mass of a binary from its inspiral phase 
to a fraction of a percent. We shall see in this section 
that the ringdown modes can determine a black hole's 
mass to a percent or tenth of a percent depending on the 
mass ratio. Therefore, the inspiral and ringdown phases 
together can shed light on how much mass is lost in the 
process of merger and how that depends on the mass ratio 
of the binary and, not probed in this study but expected 
to depend on, the spin magnitudes and orientations of 
progenitor black holes. 

When black holes merge, some of the orbital angular 
momentum goes into the final black hole. Therefore, the 
final black hole will spin in a direction different from ei- 
ther of the progenitor black holes. Independent measure- 
ments of the orbital angular momentum from the inspiral 
phase and black hole spin from the ringdown phase, could 
unravel the spin-orbit dynamics of black hole merger. 
While this is an exciting possibility, in this paper we have 
focussed only on binaries with non-spinning components. 
For such systems, the relative amplitudes of the different 
modes depend on only the mass ratio of the progenitor 
binary. For binaries with spinning black holes it is hard 
to guess how many more parameters might be required 
to characterize the relative amplitudes and hence their 
measurability. We will address this question in a forth- 
coming publication. We shall show in this section that 
one can exploit this fact to determine the mass ratio of 
the progenitor binary from the measurement of the ring- 
down mode alone. Consistency of the mass ratios from 
the inspiral and ringdown phases could offer further tests 
of general relativity. 

In this section we will explore what measurements 
might be possible by using ringdown signals alone and 
what we might learn by combining the information ob- 
tained from the inspiral phase of the signal with that 
obtained from the ringdown phase. To this end, we shall 
assume that the signals are loud and compute the co- 
variance matrix to get an estimate of the measurement 
uncertainties in the various parameters of a ringdown sig- 
nal h A (t). The covariance matrix C A m is the inverse of 
the Fisher matrix F A m given by 



F A - 

km 



dh A dh J 



d\ k ' d\ m 



X k = {M, j, q, . . .}, (39) 
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FIG. 6. Same as Fig.[5]but for LISA and the source is assumed to be at a red-shift of z = 1. Also, the mass ratio is allowed to 
vary over a larger range from 2 to 25 instead of 2 to 10. The 'steps' that can be seen around 40 x 10 6 Mq and 80 x 1O 6 M0 are 
mostly due to the LISA noise curve. 



where for any two functions g(t) and h(t) the angular 
bracket (g, h) denotes their scalar product defined by 

(g,h) = ^J °G(f)H*(f)^jr. (40) 

Here, as before, A is an index denoting the detector 
in question, Sh(f) is the one-sided noise power spectral 
density of the detector, G(f) and H(f) are the Fourier 
transforms of the time-domain functions g(t) and h(t), 
respectively, and a * denotes the complex conjugate of 
the quantity in question. The above integrals are often 
performed numerically and it is essential then to appro- 
priately choose the lower and upper limits in the integral 
so that outside this limit the integral has negligible con- 
tribution. Note that the detector noise power spectral 
density rises steeply outside a certain frequency range 
often assuring the convergence of these integrals. If a 
network of detectors is used to estimate the parameters 
then the Fisher matrix for the network is simply the sum 
of the Fisher matrices for the individual detectors: 



km 



£** 



A 

rn ' 



(41) 



A. The full parameter set 

In the case of a binary consisting of non-spinning black 
holes on a quasi-circular orbit, the effective amplitudes 
Bi m in Eq. (24 1 of the quasi-normal modes of the fi- 



where the sum is over all the detectors in the network. 



nal black hole, depend on a set of eight parameter^] 
(M, j, q, Dl, 9, ip, ip, l) : the mass M and spin magni- 
tude j of the black hole, the mass ratio q of the progen- 
itor binary, the position vector (Dl, 0, f) of the black 
hole with respect to Earth, the polarization angle ip and 
the inclination i of the black hole's spin angular momen- 
tum with respect to the line-of-sight. The phases 7£ m 
are given by Eq. ( 26 1 and they depend on the angles 
{6, tp, ip, t, cj>, (f>£ m ), where (t, 4>) are the spherical po- 
lar coordinates giving the direction in which the black 
hole quasi-normal mode is emitted in a frame fixed to 
the black hole and <fit m are the initial phases of different 
quasi-normal modes. Thus, if we consider a superposi- 



6 Recall that the final spin of the black hole is determined by the 
mass ratio of the progenitor binary and so it is not necessary to 
treat both q and j as independent. However, such a treatment 
allows us to check if the final black hole spin is consistent with 
the mass ratio as predicted by numerical relativity simulations, 
which would indeed be another test of general relativity. 
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tion of three quasi- normal modes then there will be 12 
parameters, including 022, 033 an d 4>2i- The amplitudes 
Bn m depend on eight of these parameters (exceptions are 
4>i ^22, 033 and 02 i) and phases 7^ m also depend on a 
(different) set of eight parameters (exceptions are Dl, 
M, q, and j). 



B. Measurements with a network of detectors 

Measuring all the parameters of a ringdown signal will 
require simultaneous observation of the signal in two or 
more detectors. Let us first take a look at the configura- 
tions of LISA, ET and advanced ground-based detectors. 

ET and LISA both have a triangular topology and con- 
sist of three V-shaped interferometers, with an opening 
angle of 60 degrees, rotated relative to each other by 120 
degrees. The three interferometers are completely equiv- 
alent, in terms of sensitivity, to two L-shaped interfer- 
ometers [68], with arms that are only three-quarters in 
length of the arms in the triangle. Thus, for the purpose 
of detection and measurement, we can consider ET and 
LISA to be a network of two collocated detectors. At 
least three advanced ground-based detectors (two LIGO 
detectors and Virgo) would be operating by 2015, with 
the possibility of the Japanese Large Cryogenic Gravita- 
tional Telescope joining the network soon after. Thus, 
there will be a global network of ground-based detectors 
that will be operating for a number of years from around 
2015. 

Of all the parameters, the direction to the source (6, ip) 
is the most critical and difficult to measure from the ring- 
down modes alone. However, since the ringdown modes 
we study are preceded by the inspiral phase of a binary 
coalescence, we can expect the direction to the source to 
have been measured to some degree of accuracy. In the 
case of LISA, the inspiral phase of supermassive black 
holes could last for several months to years in the detec- 
tor band. The modulation of the signal caused by LISA's 
motion relative to the source over the observation period 
will be good enough to measure the sky position (see, for 
instance, Ref. [59]). ET, together with a network of other 
detectors, advanced or 3rd generation detectors present 
at the time, should be able to triangulate the source. 
This is also true with the network of advanced detec- 
tors. Thus, we shall assume that the parameters (9, (p) 
are known, leaving 10 parameters to be measured from 
the ringdown phase. However, for very massive systems 
(depending on the detector in question), only the ring- 
down phase might be visible and for such systems it will 
not be possible to infer the location of the source with- 
out a network of detectors. Such events will not be very 
useful for testing GR. 

As expected, the relative contributions of the inspi- 
ral and ringdown phases depend on the total mass of the 
system: For lighter masses the ringdown phase makes lit- 
tle impact on parameter estimation; for heavier systems 
just the opposite is true. For systems with intermediate 



masses the contributions could be roughly equal. Such 
systems will be golden binaries with the best ability to 
test GR. 

A single detector can measure the mass M and spin j 
of the black hole by simply inverting the QNM frequen- 
cies and damping times. Additionally, each detector in 
a network would also measure three independent ampli- 
tudes B22, #21 an d -S33 and three independent phases 
722, 7^\ and 7^ — a set of 12 additional measurements 
from two detectors. Of course, the amplitudes and phases 
(as well as the time-constants and mode frequencies) are 
all expressed in terms of the 10 physical parameters and 
will not treated as independent. The counting argument 
given here shows that a set of two or more detectors al- 
lows enough measurements to fully reconstruct the ring- 
down signal. 

Therefore, one can, in principle, measure all the ten 
parameters of a QNM composed of three modes, using a 
network of two or more detectors. We have not, however, 
explored the problem in its full generality as the Fisher 
matrix that includes both intrinsic and extrinsic param- 
eters happens to be highly ill-conditioned. In such cases, 
fisher matrix is not the right approach for computing the 
errors incurred. We will, in the near future, investigate 
this problem by other means, for instance using Bayesian 
inference. For now, our goal is to see how well a subset of 
interesting parameters can be measured if, as mentioned 
earlier, we know some of the parameters from the inspiral 
phase. 

For the sake of simplicity, we shall assume that the 
phase of the different quasi-normal modes at the begin- 
ning of the ringdown are all the same and equal to zero: 
022 = 021 = 033 = 0. In this work, we have dropped them 
from further consideration so that we can focus our effort 
on the main goal of the paper, which is to show that one 
can infer the mass ratio of the progenitor binary by mea- 
suring two or more quasi-normal mode amplitudes. We 
shall, therefore, assume that the ringdown signal depends 
on the parameters (M, j, q, Z?l, t, ip, 0), seven parame- 
ters in all. A single detector can measure the mass and 
spin of the black hole from the different mode frequencies 
and damping times, as well as three amplitudes and three 
phases. Consequently, in the case of a simplified signal 
model, where we have dropped the constant phases and 
the location of the black hole from the list of parame- 
ters, we do not need a network of detectors to resolve the 
signal parameters. 



VII. UNDERSTANDING MASS LOSS, SPIN 

RE-ORIENTATION AND MASS RATIO FROM 

RINGDOWN SIGNALS 

We have computed the covariance matrix Ckm [cf. 
Eq. (39)] of ringdown signals as a function of the black 
hole mass M and the mass ratio q of the progenitor bi- 
nary. In computing the covariance matrix, we used, as 
described in the previous section, a signal model with 
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FIG. 7. Dimensionless fractional errors (Jk/Xk in various parameters as a function of the black hole mass and progenitor binary's 
mass ratio. The black hole is assumed to form at a luminosity distance of Dl = 1 Gpc and the various angles are assumed to 
be as in Table [Hi] The bottom right panel plots the correlation coefficient between the luminosity distance and the orientation 
i of the black hole's spin with respect to the line-of-sight. The general trend for the errors is to increase with increasing mass 
ratio and decreasing mass, except for Dl for which there appears to be an 'island' around q — 6. We have used the ET-B 
sensitivity curve in computing the covariance matrix. 



seven parameters: X k = {M, j, q, Dl, l, ip, (/>}. The full 
covariance matrix contains 7 variances C'kk and 21 co- 
variances C m k = Ckrm k 7^ m. The full set of results 
is too large as our covariance matrix contains 28 inde- 
pendent elements at each point in the (M, q) plane. To 
save space, we have chosen a subset of these for further 
discussion. 

We will discuss the error in the estimation <?f. — \/Cfcfe 
of 5 of the 7 parameters, X k — {At, j, q, £>l, t}, and also 



include in our discussion the covariance between the lu- 
minosity distance D^ and the inclination i of the black 
hole's spin axis with the line-of-sight. Also, it is instruc- 
tive to deal with the correlation coefficient defined as 
Ckm = Ckm/fak o~ m ), instead of the covariances them- 
selves. Correlation coefficients are bound to the range 
[— 1, 1] and capture how variation in one parameter might 
be offset by varying another. A correlation coefficient of 
for a pair of parameters indicates that they are com- 
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FIG. 8. As in Fig. [7] but for supermassive black holes observed in LISA at a red-shift of Z — 1. Note also that the mass ratio q 
in this case is allowed to vary from 2 to 25. 



pletely independent of each other and have "orthogo- 
nally" different effects on the waveform. For instance, 
the amplitude A and phase <p of a simple sinusoid func- 
tion s(t) = A sm[L>j (t — to) + f] will have a correlation 
coefficient of 0, while the phase and time offset £q are 
perfectly anti-correlated and so have a correlation coeffi- 
cient of — 1 : A change in the phase cannot be mimicked 
by a change in the amplitude but it is completely repli- 
cated by a change in the time offset. Consequently, the 
amplitude is completely independent of the phase and 
time offset, while only one of phase or time offset can be 



considered to be an independent parameteiM When co- 
variances are close to ±1, the parameters concerned will 
have large uncertainties and this, as we shall see below, 
is a major source of error for the parameterirj (Z?l, 0- 

Figures [7] and [8] plot fractional errors (i.e. a\ k /\k) 
incurred in the measurement of the parameters 
(M, q, j, £>l) and the absolute error in the parameter 
cos l. We have also plotted the correlation coefficient 



7 In fact, the Fisher matrix for the parameter set (A, ip, to) will, 
as can easily be verified, be singular. 

8 This important point was noted recently by Nissanke et al |70| . 
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cd l l, labelled in the figure as Cov(D, t). Figure [7] cor- 
responds to ET's observation of stellar and intermediate 
mass black holes in the range [10, 1000] M Q and the mass 
ratio q of the progenitor binary in the range [2, 10]. Fig- 
ure [8] corresponds to LISA's observation of supermassive 
black holes of mass M in the range [3 x 10 6 , 10 8 ] M Q and 
mass ratio q in the range [2, 25]. In the case of aLIGO, in 
most of the parameter space and for all parameters (ex- 
cept the total mass), the fractional errors are larger than 
50% and so they are not shown. Advanced LIGO will 
determine the total mass of an intermediate mass black 
hole that forms from the merger of two nearly equal mass 
black holes within 1 Gpc to within a few percents and this 
could be very interesting for some of the tests of general 
relativity to be discussed below. 

Let us recall that our results for mass ratios greater 
than 11 are based on the extrapolation of analytical fits 
to numerical simulations of binary black hole mergers 
that are only available up to a mass ratio of q = 11. In 
the next few sections we will discuss our results in the 
context of the science questions they can address. 



A. Mass ratio and component masses of the 
progenitor binary 

A key result of our study is that one can measure 
the mass ratio of a progenitor binary by observing the 
ringdown signals emitted by the black hole that forms 
from the merger. While the different mode frequencies 
and time-constants all depend only on the mass of the 
black hole and its spin magnitude, their relative ampli- 
tudes depend on the mass ratio. As we have argued be- 
fore, under certain circumstances it should be possible to 
measure the relative amplitudes of the different modes, 
thereby measure the mass ratio of the progenitor binary 
and hence deduce its component masses. 

The top right panels of Figs. [7J and [8] show the frac- 
tional accuracy with which the mass ratio can be deter- 
mined assuming that the signal is composed of three ring- 
down modes, namely (2,2), (3,3) and (2,1). ET will not 
measure the mass ratio very well in most of the parameter 
space. However, for equal mass mergers, ET should con- 
strain the mass ratio to within 5%. Remarkably, LISA 
will be able to measure the mass ratio to better than 10% 
over 60% of the parameter space studied and black holes 
that result from the merger of equal mass black holes 
enable the determination of q to better than 1%. 

Since the mass ratio of a binary is easily determined 
from the inspiral phase, its measurement also from the 
ringdown phase offers newer tests of general relativity. 



B. Mass loss to gravitational radiation 

In the process of inspiral and merger, a binary black 
hole emits a significant fraction (a few percent) of its mass 
as gravitational radiation. The total mass of a binary can 



be measured very accurately from the inspiral radiation it 
emits. Estimateaj range from a fraction of a percent (for 
an equal mass binary black hole of total mass 200 M Q at a 
distance of 1 Gpc) in the case of Einstein Telescope [7_T] to 
50 parts per million (for an equal mass binary black hole 
of total mass 2 x 10 6 M at a distance of 3 Gpc) in the 
case of LISA [26] . For most binaries observed with ET 
and LISA, the total mass before merger can be measured 
with an error that is much smaller than the fraction of 
mass that is expected to be lost in gravitational radiation 
during merger. 

What we see from the top left panels of Figs. [jjand[8] 
is that the mass of the final black hole that results from a 
merger can also be measured very accurately. ET cannot 
measure masses of stellar mass black holes very well but if 
intermediate mass black holes (and binaries composed of 
such black holes) exist, then ET will measure their masses 
to better than 1%, assuming the source is at 1 Gpc, over 
a significant range of the parameter space (q ^ 6 and 
M > 400 Mq) we explored. A black hole that results 
from the merger of two black holes each of mass roughly 
about 500 M & could be measured to an accuracy of better 
than half-a-percent. 

LISA is able to measure the mass of a supermassive 
black hole that it observes at a red-shift of z = 1 with an 
accuracy of better than 1% all over the parameter space 
that we explored. Masses of supermassive black holes 
that form from the merger of two roughly equal-mass 
(q Ss 10) black holes could be measured to an accuracy of 
0.1%. This means that from the ringdown signal alone, 
we should be able to measure masses of supermassive 
black holes even at a red-shift of z — 5 to better than 
1%. 

It would, therefore, be very interesting to compare the 
observed mass loss with the predictions of analytical and 
numerical relativity and verify if the mass loss is in ac- 
cord with their predictions. Such comparisons will put 
general relativity to new kinds of tests in the dissipative 
regime of the theory. Let us recall that the luminosity of 
a binary black hole, close to merger, could be as large as 
10 50 Js -1 , which is arguably the largest luminosity any 
physical system could have. It would be very interesting 
to test the theory when the luminosity is as large as this. 

These phenomenal accuracies with which masses can 
be measured raise the question if it is prudent to treat 
the mass of a binary to be constant in the course of its 
inspiral and merger. It might be possible to deduce the 
rate of mass loss by treating in our computation of the 



9 Note that most of the literature quotes error in the measurement 
of the chirp mass M. = M u 3 ' 5 , where M is the total mass and v 
is the symmetric mass ratio of the binary. To estimate the error 
in the total mass we have used the error propagation formula 

M > ~\M> 25 W 
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where cm v is the correlation coefficient of the parameters ln.M 
and In v. 
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waveforms the mass of the system to be a function of 
time. 



C. Exploring naked singularities 

The magnitude and orientation of the spin of the final 
black hole that results from the merger of a black hole bi- 
nary depend on a number of parameters of the progenitor 
binary: magnitudes and orientations of the spins of the 
two component black holes relative to the orbital angular 
momentum and the mass ratio of the progenitor binary. 
A spinning black hole binary has a rather large parameter 
space, six parameters more than a non-spinning system. 
Limited studies have been carried out in assessing how 
well one might be able to measure black hole spins from 
the inspiral phase of the merger of a black hole binary 
[721 173) . Numerical relativity simulations of the merger 
of spinning black holes are still in their early stages. In 
the coming years we are likely to learn a great deal about 
spin dynamics of a binary before and after merger. While 
these are important problems to be addressed in the fu- 
ture, we recall that in this paper we have only studied 
binaries comprising of initially non-spinning black holes. 

Figures [7] and [8] middle left panels, show how accu- 
rately one might be able to measure the spin magnitude 
j in ET and LISA, respectively. The accuracy here is 
not as good as in the case of the black hole mass. ET 
can deduce the final spin to within 10% over 40% of the 
parameter space and to better than 5% in 20% of the 
parameter space, for black holes that form within 1 Gpc. 
LISA, on the other hand, can measure spin magnitudes 
to better than 5% over 50% of the parameter space and 
to better than 1% over 20% of the parameter space for 
black holes that form within a red-shift of z = 1. Spins 
of black holes that form from the merger of nearly equal 
mass black holes can be measured to 0.5%. Thus, LISA 
should be able to reliably measure spin magnitudes that 
are only a few percent larger than 1. Measuring spin mag- 
nitudes to such a high accuracy will be useful in testing 
whether a merging binary results in a black hole or a 
naked singularity [17] . 



The transition from negative to positive correlation be- 
tween Dl and cos i occurs when mass ratio q ~ 6. The 
significance of this number is not clear to us at the mo- 
ment. This correlation completely destroys the accuracy 
with which the luminosity distance can be measured. For 
instance, we see that in the case of ET <td l /Dl, is in the 
range 5-10%, although the SNR in this region of the pa- 
rameter space is ~ 300. The parameter cos i in this region 
is also determined rather poorly, at about 25%. 

These numbers are better in the case of LISA: For 50% 
of the parameter space, LISA can measure the distance 
to within 10% but black holes that result from the merger 
of roughly equal mass systems might allow the luminos- 
ity distance to be determined to 1% or less. Using the 
inspiral phase alone, but with the help of sub-dominant 
signal harmonics, luminosity distance can be measured 
to a fraction of a percent in the case of LISA. Thus, it 
would be very interesting to see how similar are the dis- 
tances obtained from these different phases of the merger 
dynamics. 



VIII. PARAMETRIZATION FOR TESTING 
THE NO-HAIR THEOREM 

In this section we will consider a practical implemen- 
tation of testing general relativity using quasi normal 
modes. To test general relativity with quasi-normal 
modes, it is not necessary to consider all the physical 
parameters but only those that are necessary to fully 
characterize the shape of the signal. The ringdown sig- 
nal composed of a superposition of n quasi-normal modes 
can be written as 






cos(cj£ m t + lem), (42) 



where there are n each of the amplitudes Ae m , time- 
constants T£ m , frequencies u>t m and phases j£ m . The sig- 
nal is, therefore, characterized by a set of An parameters 
in total. 



D. Spin orientations and the luminosity distance 

The measurement accuracy of spin orientation, given 
by cos l, is shown in the bottom left panels of Figs. [7] and 
|8j cos i can be measured to within 10% in about one- 
third of the parameter space in the case of ET and to 
better than 5% in about one-half of the parameter space 
in the case of LISA. The dynamics of spins before and 
after merger could be relevant in understanding the x- 
shaped radio galaxies [74"] . 

The spin orientation of a black hole is very strongly 
correlated with the luminosity distance. We see from 
the bottom right panels that the correlation coefficient is 
close to either +1 or —1 in most of the parameter space. 



A. Maximal set 

The most exhaustive test of the no- hair theorem would 
be to treat all An parameters to be independent. Of the 
An parameters, only the 2n mode frequencies and time- 
constants would facilitate the test, the others should be 
retained in order to fully capture the covariances and 
variances in the 2n test parameters. The consistency 
among every mode frequency and time constant makes 
the test more stringent but the presence of a large num- 
ber of parameters (when n is greater than 2) weakens the 
test. The reason for the latter is that a model with too 
many parameters will/should be penalized for its flexi- 
bility by any carefully formulated test. In a covariance 
matrix formulation of the test, this will be reflected by 
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FIG. 9. Curves of constant mode frequencies and time-constants in the (M, j)-plane obtained with Eqs. ( 43 1 for the (2,2) mode 
and similar equations for other modes. Top panels correspond to the case when all measured values are exactly as predicted 
by black hole perturbation theory. The bottom panels correspond to the case where the (3,3) and (2,1) mode frequencies and 
time constants differ from the GR value by 10% but the (2,2) frequency and time-constant are as in GR. An interesting thing 
to note is that for a large range of values of the black hole spin, the curves of constant T22 and T33 almost overlap, providing 
the opportunity for a more accurate test of the no-hair theorem, (see the text for details.) 



large variances in T£ m and ui m (which are our test pa- 
rameters) and in a Bayesian model selection a model with 
a larger number of parameters will suffer from having a 
large evidence. 

How would the test work in practice? For each mea- 
sured test parameter (and the associated error in its mea- 
surement), one could draw a curve (or a band including 
the error) in the (M, j) plane, by using their expressions 
in general relativity in terms of the mass and spin of the 
black hole. If the curves/bands fail to intersect at a single 
point/region in the (M, j) plane then that would invali- 
date general relativity or, alternatively, indicate that the 
object is not a black hole. The maximal set could be 
weak due to the large variances of the various param- 
eters but strong because many different bands have to 
pass through the same region. 

Let us illustrate how the test works with some exam- 
ples. Let us suppose our signal model consists of a su- 
perposition of {£, to) = (2, 2), (2, 1), (3, 3), (4, 4) modes. 
In this case the maximal set contains 16 parameters, of 
which 8 mode frequencies and time-constants are the test 
parameters. Let us denote by &i m and fi m , the values 
of the mode frequencies and time-constants measured by 



projecting the data onto a superposition of quasi-normal 
modes as in Eq. (42). For each measured parameter, we 



can construct an equation in (M, j) using the relation be- 
tween the parameter and (M, j) given in Eqs. (27)-(30): 



W22 = ^ [1-5251 



T22 



2 

w 22 



1.1568(1 -j) ' 1292 ], 
[0.7000 + 1.4187(1 _j)-°-4990] 



(43) 

and similar equations for other modes. (Note that we 
have used /22 = M0J22 and 2 Qt m = Ti m U)t m in rewriting 
these equations as dimensionful quantities.) Measure- 
ment errors can be folded into the analysis by using 



W22 i a U22 
T22 ± o> 22 = 



= A, [1.5251- 1.1568(1 -j) ' 1292 ], 



2 

W22 



[0.7000 + 1.4187(1 - j 



,-0.49901 



(44) 

with similar equations for other modes. 

For a 5 x 10 6 M black hole, the mode frequencies 
and time-constants for three different spin values j = 
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0.1, 0.7, 0.9, 0.99, are given in TablefEVl If the ringdown 



TABLE IV. The frequencies Fg m — <^ m /(27r) (in mHz) and 
time-constants T£ m (in s) of the first four dominant modes for 
a 5 x 10 6 Mq black hole of different spin magnitudes j. 



3 


F 22 


F 21 


F 33 


F 44 


T22 


T21 


T33 


T44 


0.10 


2.48 


2.43 


3.98 


5.36 


282 


277 


269 


266 


0.70 


3.46 


2.96 


5.48 


7.44 


303 


301 


295 


290 


0.90 


4.30 


3.30 


6.70 


9.06 


383 


355 


379 


375 



0.99 5.73 3.66 8.60 11.4 823 559 828 837 

signal is consistent with the formation of a black hole 
with spin magnitude, say, j — 0.7 (j = 0.99) then the 
mode frequencies and time- const ants would be precisely 
as in the 2nd (respectively, 4th) row of Table IV mod- 
ulo the errors in their measurement. Therefore, curves 
defined by uj em (M,j) = uj lm and n m (M,j) = T im will 
all meet at a single point in the (M, j) plane as in the 
upper two panels of Fig. |9l the point of intersection giv- 
ing the mass and spin magnitude of the black hole. The 
left panel corresponds to the formation of a black hole of 
spin magnitude j — 0.7 and the right panel to j = 0.99, 
in both cases M = 5 X 10 6 M e . If, however, one of the 
mode frequencies, say W33 is different from the general 
relativistic value by 10%, then the corresponding curves 
would fail to meet as shown in the lower two panels of 
Fig. [9] 

If general relativity is true then some of the curves lie 
almost one on top of the other (e.g. w 2 2 = const, is almost 
identical to W21 == const, in the j — 0.7 case and T22 = 
const, is identical to T33 = const, both in the j — 0.7 and 
j = 0.99 cases) but even a slight departure from general 
relativity will lead to big departures as demonstrated by 
the lower panels of the same figure. In any measurement 
the parameters are subject to statistical and systematic 
errors that must be folded into the analysis which will be 
taken up in a forthcoming study. 

Do we need to treat all time-constants and mode fre- 
quencies to be independent in a test of the no-hair the- 
orem? We shall argue below that it is not necessary to 
treat all 4n parameters to be independent; in fact, we 
shall see that the parametrization is not unique, offering 
a lot of flexibility in testing GR. 



B. Minimal set 

The minimal, or the simplest, model would consist of 
the smallest number of parameters needed to check the 
consistency between the modes as predicted by the no- 
hair theorem, yet large enough to capture all the vari- 
ances and covariances between the parameters of inter- 
est. Since the mode frequencies and time-constants are 
all determined in GR by the mass M and spin magnitude 
j of the black hole, the smallest number of parameters 
required to test GR would be three: these could be any 



three time-constants or two time-constants and one mode 
frequency, and so on. It would then be necessary to ex- 
press the other mode frequencies and time-constants in 
terms of any two of the three parameters that were taken 
to be independent. Two of the three independent param- 
eters could be used to solve for (M, j). One could then see 
if the measured value of the third parameter is consistent 
with its predicted value based on the values of M and j. 
Let us note that without a prior knowledge of the ampli- 
tudes Ai m and phases 7^ m it will not be possible to mea- 
sure the chosen time-constants and mode frequencies as 
they would induce covariances that cannot be neglected 
in estimating the errors incurred in their measurement. 
Thus, we (tentatively) conclude that the minimal set re- 
quired for a signal model with n modes would be In + 3. 
However, the three test parameters can be chosen in any 
way one wishes but choosing only three assures that the 
error in their measurement is the smallest. In this sense 
the minimal set could be a very stringent test of GR. 

In reality, of course, the amplitudes of different 
modes are determined by the physical parameters 
(M, j, q, Dl, 9, (p, tp, l) (cf. the discussion at the be- 
ginning of this section). We can assume the parameters 
(<7, Dl, 8, f) to be known from the inspiral phase and 
(M, j) to be determined by the time-constants and/or 
mode frequencies of the ringdown signal. This leaves the 
two angles (tf>, 1). Thus, when n is greater than 2, it is 
not necessary to consider all mode amplitudes to be in- 
dependent but just two of them. Thus the minimal set of 
parameters to be considered for testing the no-hair the- 
orem is n + 5 (n phases, 7£ m , the 3 test parameters and 
the two angles tp and t). 



IX. CONCLUSIONS AND FUTURE WORK 

In this paper, we have explored what information can 
be extracted from a black hole's ringdown signal, wherein 
the perturbation is caused by the tidal deformation pro- 
duced during the merger of two non-spinning black holes. 
To this end we used numerical simulations of the late in- 
spiral and ringdown to estimate the relative amplitudes 
of the various modes excited. The simulations consisted 
of initially non-spinning black holes in quasi-circular or- 
bits for several mass ratios, ranging from 1:1 to 1:11. 

We find that several modes have large enough lumi- 
nosity — or signal-to-noise ratio — to be detectable in 
LISA and ET. Specifically, in order of decreasing power, 
modes (2,2), (3,3), (2,1), (3,2), (4,4), (5,5), (4,3), (6,6), 
(5,4) and (4,2) for LISA and the first four to five modes 
for ET, have significant luminosities. Note that not all 
of these can necessarily be resolved, but it probably is 
the case for the first three (see below). In the analysis 
though, we decided to include only (2,2), (3,3) and (2,1), 
mainly because the available data were most accurate for 
those modes. 

We argued that the ringdown signal depends on the 
mass ratio of the progenitor binary and that this can be 
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measured, with an error that is estimated from a Fisher 
matrix analysis. We showed how the luminosities change 
with the mass ratio. Indeed, by constructing fits to the 
mode amplitudes in terms of the mass ratio, we were able 
to include this effect in the analysis and estimate the er- 
rors involved in the measurement of various parameters. 
An important issue was to determine the epoch when 
the ringdown phase starts, so as to evaluate the rela- 
tive amplitudes at that point. This was taken to be the 
point where each mode's frequency stops having an up- 
ward trend. The epoch at which the peak luminosity is 
reached is slightly different for different modes, but the 
mode amplitudes were all measured at a time 10M after 
the peak luminosity of the 22 mode. 

We computed the measurement errors of a number of 
other parameters for ringdowns observed in LISA and 
ET. These include the black hole mass, its spin, lumi- 
nosity distance and inclination angle and how these vary 
with the final BH mass and the binary mass ratio. If 
we do not consider initial BH spins, the inclination angle 
can be measured solely by observing the three most dom- 
inant modes of the ringdown waveform, assuming that 
these can be resolved. 

Together, LISA and ET will be able to provide ample 
evidence for the distribution of supermassive, intermedi- 
ate mass and stellar mass black holes, for a large part of 
the known Universe. For most of the parameter space, 
the reach of LISA for ringdowns is z ~ 6, while for ET 
at least z ~ 0.8. Also, by being able to measure the 
mass ratio, hints on the merger history and formation of 
black holes of a large range of masses could be inferred 
by studying ringdown signals. 

A practical implementation for testing the no-hair the- 
orem and deciding the nature of the compact object that 
results from the merger, was presented, illustrating sev- 
eral key components. We started by providing a general 
framework, which is based on the number of parame- 
ters necessary to apply the test. Specifically, using n 
modes, a generic test will use the 2n frequencies and time- 
constants. This test is implemented by plotting the cui m 
and T£ m curves on the mass-spin plane of the final black 
hole, where all the curves should intersect inside the same 
region if the object is a black hole. A key point was that 
some of the curves are special, as they almost overlap 



and can be thus used to check for small deviations from 
general relativity. 

Future work should extend the study to include nu- 
merical simulations of initially spinning black holes, as 
this corresponds to a more realistic scenario. Addition- 
ally, there is effort to produce more accurate numerical 
simulations so that less dominant modes could be studied 
when the mass ratio is large. 

The current study has not investigated the question 
of decomposition of modes of a ringdown signal in real 
data. The SNRs of different modes, especially in the case 
of LISA, suggest that it should be possible to resolve the 
modes and carry out the proposed tests of general relativ- 
ity. However, a more in-depth investigation needs to be 
done, for instance by using a Bayesian model selection to 
discriminate between different models. Given some prior 
information, the relative probability of two different mul- 
timode ringdown waveforms - injected in Gaussian and 
stationary noisy data - could be computed. 

Finally, a Fisher matrix-based parameter estimation 
approach (such as the one presented in this work) may 
not be robust or accurate enough, especially when the 
parameter space is large and the signal-to-noise ratios 
are low [75]. For this reason, future studies should also 
involve parameter estimation in the context of the afore 
mentioned Bayesian analysis, as well as a Bayesian ap- 
proach to a test of the no-hair theorem. In the latter case 
the posterior joint probabilities of two different models, 
given the 'initial data', could be compared. The first 
hypothesis would be that the observed object is a black 
hole, while the second that either GR is incorrect or the 
merged object is not a black hole. 
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